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Abstract 

The asymptotics of the generic second Painleve transcendent u{x) as 
\x\ — » co, argrr. £ (f^ f(' c + !))> ft = 0, 1, . . . , 5 is found and justified via 
the direct asymptotic analysis of the associated Riemann-Hilbert prob- 
lem based on the Deift-Zhou nonlinear steepest descent method. The 
asymptotics is proved of the Boutroux type, i.e. it is expressed in terms of 
the elliptic functions. Kapaev-Novokshenov's explicit connection formu- 
lae between the asymptotic phases in the different sectors are obtained as 
well. 



1 Introduction 

The large argument asymptotic behavior of the Painleve functions of the first 
and second kinds were first described in the classical work of P. Boutroux Q . 
In particular, Boutroux found a general 2-parameter elliptic (i.e. expressed in 
terms of elliptic functions) asymptotic formulae for the Painleve functions. 

In 1988, using the multiscale expansion method, N. Joshi and M. Kruskal 
p6[ revised the result of Boutroux and found the phase shift in the Boutroux 
elliptic ansatz. Specifically, Joshi and Kruskal showed that in the case of the 
second Painleve equation, 

d 2 u o 

— = 2u 3 + xu, (1.1) 

ax* 

the large \x\ asymptotics of its generic solution u(x) is given by the following 
equation, 

u(x) ~ civ^sn (-C2X 3 / 2 + fJk)\ + uui2\D). (1-2) 
o 

Here the periods u)j, j = 1,2 of the elliptic sine, the multipliers C1.2, and the 
module D are (transcendental) functions of argx, while the parameters n, v are 
constant in the interior of any sector argx £ (-|fc, -|(fc + 1)) , k = 0, 1, . . . , 5. 



The connection formulae between the constants /i, v in the ansatz (1.2) 
corresponding to the different sectors were found by V. Novokshenov J2(| and 
A. Kapaev pO, Ol| with the help of the isomonodromy deformation method 



which was introduced in the very beginning of the 80-s by H. Flaschka and 



A. Newell § and by M. Jimbo, T. Miwa and K. Ueno [17|, gj| |L9|. 

The isomonodromy deformation method associates with a given Painleve 
equation an auxiliary system of linear ODEs with rational coefficients whose 
monodromy data are the first integrals of the Painleve equation. Hence the 
asymptotic evaluation of the solutions of the Painleve equations reduces to the 
asymptotic solution either of the direct or of the inverse monodromy problems 
for the auxiliary linear system. The relevant asymptotic scheme based on the 
asymptotic analysis of the direct monodromy problem was suggested in [jlJJ 
and was used there and in the subsequent works of several other authors for 
the evaluation of the connection formulae for the Painleve transcendent. We 
refer the reader to review Jl2| for more on the history and apparatus of the 
isomonodromy method. 

The approach of J|(| and [^o| . is based on an extension of the direct 
monodromy problem technique of |14| to the asymptotic analysis of the Painleve 
functions in the complex domain. The further development of this scheme was 
made in the works of A. Kitaev and A. Kapaev and A. Kitaev [^3| which 
were also devoted to the elliptic asymptotics of the Painleve transcendents. 

In the beginning of the 90s, P. Deift and X. Zhou proposed a nonlinear 
version of the classical steepest descent method for oscillatory Riemann-Hilbert 
problems (^] . This allowed to develop an alternative to |l4| asymptotic approach 
for solving the Painleve equations based on the asymptotic solution of the rele- 
vant inverse monodromy problems, which can be formulated as the oscillatory 
matrix Riemann-Hilbert (RH) problems. This approach has the advantage of 
not using any prior information of the solutions of the Painleve equations, and 
it has already been applied to the trigonometric and exponentially decreasing 
asymptotics of the Painleve functions on the real axes 0, 13 1. 



The goal of this paper is to extend the Deift-Zhou method, in the case of 



the second Painleve equation (IT), to the complex domain. As a result we will 
obtain, rigorously and without any prior assumption crucial for and [ po[ 
the Boutroux type elliptic asymptotics together with the explicit formulae for 
the relevant phase shifts directly from the associated Riemann-Hilbert problem. 
The plan of the paper is as follows. In Section 0, we remind the setting 



the RH problem related to equation (1.1) and formulate our main theorem 
(Theorem 2.1). In Section [3], we transform the RH problem to the one posed on 
the graph consisting of the anti-Stokes lines of a certain function g(z) emanating 
from its critical points and connected by three finite Stokes lines. In Section |], 
we explicitly construct the function g(z) mentionec^]. In Section^, we construct 
the solution of the reduced RH problem on the finite Stokes lines of g(z) in 
terms of the Jacobi theta functions. In Section ^, we construct the approximate 
solutions of the RH problem near the critical points of g(z) in terms of the Airy 
functions. In Section |?], we collect all the functions together and prove that 

1 In the most recentapplications of the nonlinear steepest descent method to the orthogonal 
polynomials (see e.g. |4l) this step corresponds to the construction of the so called equilibrium 
measure 
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this combined function approximates the solution of the original RH problem. 
In Section |^, we extract from the approximate solution of the RH problem the 
approximate solution of the second Painleve equation. 

The authors dedicate this paper to Barry McCoy on the occasion of his 
60th birthday. We are especially happy to do this since it was a remarkable 
1977 paper of McCoy Tracy, and Wu |^5j where it was demonstrated for the 
very first time that an explicit and complete connection formulae for a Painleve 
transcendent are possible. 

2 RH parametrization of the Painleve functions 
of the second kind. Formulation of the main 
theorem 



According to Flaschka and Newell || (see also Q and [jl4 



) the Riemann-Hilbert 



problem associated with the second Painleve equation (1.1) consists of finding 
the piecewise holomorphic 2x2 matrix- valued function ^f(X) of the complex 
variable A such that 

1) \I/(A) is holomorphic for A G C \ U{7fc}, where -fk are the rays jk — {A G 
C : arg A = -| + -|(fc — 1)}, k = 1, . . . , 6, oriented from zero to infinity, and 

■$(X)e e(X ' x)(73 -> /, A^oo, 0(X,x) =i(^-\ 3 + x\), (2.1) 
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2) on the rays "fk the jump conditions hold 

*+(A) = *_(A)S fe , Ae 7fc , (2.2) 

where 

s — UL :)■ M» (2 - 3) 

and the parameters Sk do not depend neither on x nor on A and satisfy 
the constraints 

Sfc+3 = -Sk, si - s 2 + s 3 + S1S2S3 = 0. (2.4) 

Here, ^+(A) and 5'-(A) are the respective limits of ^(A) to the left and 
to the right of 7. 

The RH problem is depicted in Figure 1. 

It is worth noticing that the first of the relations ( |2.4| ) and the central sym- 
metry of the RH problem graph imply the symmetry equation 

*(A) = <7 2 *(-A)ct 2 , <r 2 =(°. ~*Y (2.5) 
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Figure 1: The RH problem graph 7 for the ^-function. 



Also, the special triangular structure of the matrices ensures that the product 
e~ B<J3 Ske e ' J:i — » I as A — > 00, i.e. the well-posedness of the RH problem. 

The piecewise holomorphic function ^(A) is in fact a collection of the canon- 
ical fundamental solutions of the linear 2x2 matrix ODE 



— = {-i(4A 2 + x + 2u 2 )a 3 - 4u\a 2 - 2vai}V, 
cLX 



(2.6) 



a x = 



1 

1 



where v and u are (meromorphic) functions of x such that v — u x and u(x) 
satisfies the second Painleve equation (1.1). The jump matrices Sfc of the RH 
problem have the meaning of the Stokes matrices for ( |2.6| ) and, b y co nstruction, 
do not depend on x. Therefore, the secon d P ainleve equation (1.1) describes 
the isomonodromy deformations of system (2.6). 

The Painleve function u(x) admits the following representation in terms of 
the solution 'J (A) of the RH problem, 



u[x) = 2 lim ( A*i 2 (A)e 

A — M 



-B(\,x) 



(2.7) 
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Moreover, the respective map, 



{s} = {s k : s k+3 = -s k , si - s 2 + s 3 + sis 2 s 3 = 0} 1 > {solutions of (|l.l|)}, 
is a bijection. Hence we adopt the notation 

u(ar) = u(x|s) 



for the second Painleve transcendent. 

The following important symmetry properties take place, 



u{x\s) = e^u^xe^s 21 ), 



(2.8) 
(2.9) 



u(x\s) = u(x\s*), 

where s 1— > s , e? G Z and s 1— > s* denote the automorphisms of the set of 
monodromy data {s} defined by the relations 



(s d )j := s j+d , 



and 



(s*)j : = 84-3- 

Symmetries (B]q) and ( |2.9| ) show that to obtain a complete asymptotic descrip- 
tion of the function u(x\s) as x — > 00, x £ C , it is enough to find, in terms of 
s, the asymptotic behavior just in one of the sectors (fk, + 1)) , say in the 
sector (4p 7r) . The asymptotics of u(x\s) in this sector is then of our principal 
interest. 

Let us introduce a complex valued function h — >c(ip), < ip < ^, which 
satisfies the system of transcendent equations 



(-|^) 3/2 (-2^ + + 

p iV s 3/2 

' ((^-ljjf + Ci + **)£;) 



Re 
Im , 

.1 + x 2 

where K and E are the standard complete elliptic integrals 



= 0, 

= 0, 



(2.10) 



f 1 dz 
K = K(x) = / , K' = K(x>), 

Jo V(l-^)(1-^ 2 ^ 2 ) 



E = E{x) = J 1 ]/ 1 1 ff dz, E' = E{^\ 

x = \J\ - >f 2 , < argx:' < ^, < arg(l + x 2 ) < 0. 

One can show that system ( [2.10 ) has a unique smooth solution normalized by 
the conditions. 



2vr 



7T =>• X — * 0. 
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-- < argx< 0. 

The main result concerning the asymptotics of u(x\s) for complex x can be now 
formulated as the following theorem. 



Theorem 2.1 Let 

Then 



S 3 ^0, 1-S1S3^0. (2.11) 



1/2 (2i x 3/2 2iK 
VI + >r \o Vl + >r TT 

+ ^ln(l_ SlS3 ) + 0(x- 3 / 2 )|>f), (2.12) 
2n 

z->oc, argx = (y3 + 0(x" 3/2 ), ¥>e(— ,7r), (2.13) 

o 

and t/ie asymptotic is uniform on the cheese-type domain 

{xeC: |arga;-^| < C\x\~ 3/2 } n D £ , 

where D e is the complement to the union of the e -neighborhoods of all the poles 
of the indicated elliptic function. 

This theorem is due to V. Novokshenov j2(J and A. Kapaev ^0|. Formulae 
( pTC ) and (2.12) (without the indicated expression for the phase) constitute the 



classical Boutroux ansatz. In this paper we present a new proof of this theorem 
based on the Deift-Zhou nonlinear steepest descent method. This proof is free 
from any prior assumptions, which were crucial for p6| and EG] , including the 



assumption of the solvability of the RH problem under the conditions (2.11) 
The latter we will obtain en route. 



We complete this section by noticing that Theorem 2.1 together with the 



symmetry relations ( |2.8| ) and (2.9) yield the following asymptotic description 
for the rest of the sectors. 

Corollary 2.1 // S2+21 7^ and 1 + S2+21S3+21 7^ for I = 0, ±1, then the 
leading asymptotic term of the Painleve function in the interior of the sector 
argx e (— ^-l; | - 4p) is elliptic. 

If S2+21 7^ and 1 + s\ + 2iS2+2i 7^ for I — 0, ±1, then the leading asymp- 
totic term of the Painleve function in the interior of the sector argx G (— ^ — 
^-l; —jrl) is elliptic. 

The phase shift of the elliptic asymptotics in each of the sectors is given by 
the expression, 

2iK K< 

in sk ± — m sk 1 , 

where the choice of the monodromy parameters sk o,nd sk> is indicated in Fig- 



ure 2. Together with the monodromy equation (2.4) this constitutes the connec 



tion formulae for the Painleve function u{x\s) in the complex domain. 



G 




Figure 2: Combinations of the Stokes multipliers whose non-triviality yields the 
elliptic asymptotic behavior of Painleve function in the corresponding sector 



Remark 2.1 A complete description of the asymptotic behavior ofu(x\s) which 
includes the "Stokes rays" arga; = ^k is given in In fact, the complete 

description of the asymptotic behavior of the solutions of the second Painleve 
equation obtained in JH|/ covers the general case of the latter, i.e. the equation 

d 2 u 

— r = 2u 3 + xu + a, a£C. (2.14) 

aar 

3 The transformation of the RH problem 

Let us scale the independent variable A according to the equation 

A = {e- itp x) 1/2 z, (e- tv x) 3/2 = t, <p = const e tt), (3.1) 
Ret — > +oo, |Imt| < const. 

The stationary points of the phase 9 — it(^z 3 + e Iv z) are independent of 
t, and formally, the problem can be treated in the very same way as the RH 
problem which is dealt with in ref. ^ |l5| where the case x — ► — oo, i.e. ip = n, 
was considered. Indeed, applying literally the same arguments as in ref. (Is] , one 
can easily replace rays jk by the anti-Stokes lines 7^ = {lm9 = 0} indicated in 
Figure 3. Points z± are the stationary phase points, z± — =F|e lv> / 2 , and, unlike 
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Figure 3: The anti-Stokes lines for the phase 9(z) — it(^z 3 + e lv z). 



the case ip = n, the exponent 9{z) has nonzero real parts at z = z±: 

Re9(z±) — > ±^ cos i — > +00; 

since </? € (|7r,7r), for sufficiently large i, 

Re<9(z+)<0, Rc(9(z_)>0. (3.2) 

Unlike the case ip = ir, now, after the normalization of the RH problem (i.e. 
"J" — ► y = \E'e 6 ' < ' z * |cr3 ), only half of the jump matrices, namely the matrices corre- 
sponding to the lines 7+ , 7^ , 7^ , 7^" , approach the identity as t — > 00 uniformly 
for z e 7^. The other half of the jump matrices, i.e. the matrices associated 
with the lines 7^ , 73" , 7^ , 7^ , "explode" as t — > 00 in the finite neighborhoods 
of the points z±. Therefore, we have no reason to expect that the RH problem 
under consideration can be reduced to the Weber-Hermite model RH problems 
associated to the stationary phase points z±. Thus, the principal idea of the 
asymptotic analysis in the case ip = ir fails as ip £ (^ L ,7r). 

A regular procedure that allows us to overcome the indicated above obstacles 
is based on the following construction. Firstly, we "split" the stationary phase 
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points z± into four new points 



zi, z 2 = -Zx, z a , Zi = -za, 

where z±, z 3 lie in a neighborhood of the stationary phase point z+, while Z2, 
Z4 are associated in the same way with the stationary point z_. Secondly 
we assume that there exists an analytical function g(z), which satishes the 
conditions: 



a).g(z) = ? ;(|z 3 + e^z)+0(i), 



as z 



b) points ±zi, ±z 3 are the branch points of g'(z) of order 2 and, in their 
neighborhoods, 

g'(z)=cl 3 (z±z li3 ) 1 / 2 + ... , 

with cf 3 7^ 0. 

c) Reg(±zi y3 ) = 0. 

d) 9 7 0onC\{±z 1|3 }. 

In fact, conditions a)-d) constitute a certain problem of the theory of an- 
alytic functions. Simultaneously, they are the restrictions on the points ±£1,3 
themselves. In Section [|, we will solve the problem by the use of an elliptic 
integral associated with the elliptic curve w 2 = (z 2 — z 2 )(z 2 — zf). We see then 
that condition c) is just the classical Boutroux equations, which always appear 
in the asymptotic analysis of the Painleve equations. It also should be empha- 
sizing that a flexibility we gain introducing four parameters, i.e. points ±zi,3, 
plays the crucial role for solvability of the problem a)-d). 

Having assumed the existence of the function g(z), our main idea is to re- 
place the original phase 8 = ii(fz 3 + e lv z) by the function g(z). Because of 
the asymptotic equation a), this replacement preserves the normalization con- 



dition (2.1). Moreover, due to condition c), we eliminate the main trouble, i.e. 



inequalities (3.2). Nevertheless, there is a price: now, we have to deform the 
initial RH problem to the RH problem formulated for the anti-Stokes lines 7^) 
corresponding to the function g(z): 

lmg(z) = Img(zk). (3.3) 



The anti-Stokes graph defined in ( |3.3| ) is more complicated than the 0-graph 
depicted in Figure 3. Properties a), b) and d) lead to the g-anti-Stokes lines 

~(k) 

7s as shown in Figure 4. Our main task for this section, is the transformation 
of the initial RH graph presented in Figure 1 into the one given in Figure 4. 

The notation 7! for the anti-Stokes line used here refers to its starting point 
Zk and to the anti-Stokes ray j s as being asymptotic for the line considered. 

As it is easy to see, the transformation problem is not trivial and can not be 
reduced to the simple bending of the graph branches due to the very different 
topological properties of the graphs presented in Figure 1 and Figure 4. In 
particular, the RH problem graph 7 for the ^ function is connected, while the 
anti-Stokes graph 7 is disconnected. The most we can do with the original 
problem by means of trivial bending of the branches is shown in Figure 5. 
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Note that after such deformation, a new boundary curve between the third 
and the sixth domains arises where the jump is described by the ratio of the 
sixth and the third canonical solutions: 

* 6 (*) = f3(z)S 3 S4S 5 , (3.4) 

where 

c c c ( 1 + SlS 2 -.si 

->3->4->5 — i 

V - s i I-S1S3 

Let us consider the function S(z) depending on Zk, k — 1,3, and some 
parameters V[, I = 1, 2, 3, 

which is defined on the complex plane cut along the broken line [23,24] U 
[zi, —Zi] U [—21, —23]. The function is fixed by condition 

S(z) — > 1, (3.6) 

z— >oc 
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Figure 5: The continuously deformed RH problem graph. 



and has the jumps on the cuts 

S-(z) = 5+(z)e 2 ™\ ze(z 3 , Zl ), 

6_(z) = S + (z)e 2 ™*, zeizu-zj, 

5-{z) = S + (z)e 2 ^ , z e(- Zl ,-z 3 ), 



(3.7) 



as it is shown in Figure 6. 

The next step is the most important for the transformation of the original 
RH problem. Let us introduce the new function 



(3.8) 



This new function does not differ from the original one at infinity, but its jumps 
are described by the new matrices: 



(3.9) 



The new RH graph is shown in Figure 7, and the new jump matrices S are given 
by the equations, 
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Figure 6: The graph of the scalar RH problem solved by the function 8(z). 



S3S4S5 — 



(1 + s lS2 )e 2 ™ 2 -sidle- 2 ™ 2 
- Sl 5Z 2 e 27:W2 (1 - sis^e' 2 ™ 2 



Now, we are going to explain what this change is made for. The point is 
that some of the new jump matrices can be factorized in such a way, that the 
simple transformation of the RH graph allows us to approach the anti-Stokes 
lines indicated in Figure 4 (cf. Indeed, let us assume that the generic 

conditions hold true, 

s 3 ^0, 1- Sl s 3 ^0. (3.11) 

If we then put 

e- 2mU3 = s 3 , (3.12) 
then it is easy to check the equation holds: 



S 3 = S^SaS4» (3-13) 



where 



s-(s *-?')■ -*-)•*-(; ^').(-) 

Similarly, if we put 

e 2 " Wl =s 3 , or V! = -V3, (3.15) 
then the factorization arises 

S^ = sfs B S5, (3.16) 

with the matrices 

Si = ((^ s V i 1)' Sb= (-^ 2 5 ~) ' S5 = ((^V i 1 

(3.17) 
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Figure 7: The RH problem for 



At the same time, if we put 

e 2 ™ 2 = l-s lS3 , (3.18) 

tllCn .U.L , 

S3S4S5 = S S (3.19) 

with 

The "intermediate" RH graph is shown in Figure 8, and the final graph is 
presented in Figure 9. 

Here, the jump matrices are 
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Figure 9: The final RH problem graph for "J. 



8-(z) = 5 + (z)(l-8is 3 ), z£{z x ,-zi), (3.26) 
S-(z) = 6+(z) ■ ±, ze(-z 1 ,-z 3 ). 

In this form, the RH problem graph for the ^-function coincides with the 
anti-Stokes graph 7, see Figure 4, except for the additional segments [2:3,21] 
and [-21,-23]. 

Now, we do not need t he function 6 anymore and can return to the original 
^-function by the use of (3.8), 

tf = *<T 3 . (3.27) 

The calculations above show that the solution of the second Painleve equation 
can be obtained from the ^-function with the jumps on the graph shown on 
Figure 10. The new RH problem, which replaces the problem ( |2.l| )-(fO|), can 
be formulated as follows. 

1) Find a piecewise holomorphic matrix- valued 2x2 function of the complex 
variable 2 with an irregular singular point at infinity where it has the 
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Figure 10: The ultimate RH problem graph for 



exponential asymptotic behavior of the type ( |2.l| ) 

= (7 + 0(z- 1 ))e^ te(2)<T3 , 0(z) = i(^ 3 + e i¥! z), z -» oc; (3.28) 

2) on the piecewise smooth oriented graph 7, see Figure 10, the jump condition 
holds true, 

*+(z) = *_(z)S, Ae 7 . (3.29) 

Here, *+(z) and ^f-(z) are the limits of the function *B(z) on the contour 
7 to the left and to the right, respectively, and the piecewise constant 
matrix S is described by the following equations, 

s?' = (i :)• s„-' = (; ■■), s?' = (i (mo) 

(3.31) 
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s (2) = /i 
4 1 1 



, S 



(2) 



1-sisg 



s„ = 



33 , s 3 



1 o 

S3 1 



_j_ 

S3 



;(4) 



1 



(3.32) 
(3.33) 



Sn = l-«i»3 



I-S1S3 



The parameters s/. satisfy the constraint (2.4), 

si - s 2 + S3 + sis 2 s 3 = 0, 
and do not depend neither on t nor on z. 



(3.34) 
(3.35) 



Remark 3.1 Possibly, given the jump contour 7, the resulting jump conditions 
(3.3L )-(3.34) can be reconstructed directly from the jump conditions (2.3) for 
the contour 7, see Figure 1, using, instead of the "ghost" function 8{z), the 
system of algebraic relations for the matrices and 



4 Construction of the function g(z) 

In this section, we construct the function g(z) described in Section 3. Consider 
the function 

w(z) = yj(z*-z?)(z*-z£), 4 + 4 = -\e* p , (4- 1 ) 

defined on the Riemann surface T, see Figure 11, pasted of two complex planes 
cut along the intervals [23,21] and [—21,-23], while the branch of the root on 
the upper sheet is chosen in such a way that 

w = z 2 - -(21 + 2§) + 0(z~ 2 ) as z^oo + . 

Here oo + denotes the infinite point of the upper sheet. We choose the basis 
{a, b} of the group Hi(T) as it is indicated in Figure 11. Introduce the Abelian 
(elliptic, in fact) integral of the second kind 

g(z) =U [ w{z)dz. (4.2) 
Jo + 

We shall consider g(z) as a single valued function defined on the complex plane 
C * cut along the sum of the intervals [23, z\] U [21, —21] U [—21, —23], moreover 
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Figure 11: The Riemann surface of the function w(z). a and b mean the basic 
cycles of the surface. 



the starting point z = 0+ of the integration path lies on the left side of the 
oriented cut [z\, —zi]. Introduce the periods of (4.2), 



A = AiJ a = Aif w(z) dz , B = AiJ b = Ai f w(z) dz . (4.3) 

Ja Jb 

Taking into account the symmetry w(—z) = w(z), we find that 

g{z) + g(-z) = Ai( f +f +/+/ +/ )w(z)dz = A. 
Thus the function 

9(z) = 9(z) - \A, (4.4) 

considered on C*, is odd, 

gi-z) = -g{z)- 



Using definition (14), we see that properties b) and d) of g(z) are satisfied. 
To find the asymptotics of g{z) as z — > oo, introduce the regularized integral 

AzJ* (w(z) - z 2 + \{z\ + zi)) dz = g(z) ^(^ 3 + e*z). 
On the other hand, again using the symmetry w(z) — w{—z), it is equal to 
2i( £~ + [)H Z ) - z2 + \( z l + 4)) dz = 
= 2i (f(w(z) -z 2 + hz\ + zi)) dz + 0{z- r ) = \a + Oiz- 1 ). 

Thus, 

g(z) = i(^z 3 + e ztp z) +0(z~ 1 ) as z -» oo, (4.5) 

and we have established the property a) of the function g(z). The property c) of 
the function g(z), i.e. He g(±z 1,3) = 0, is equivalent to the system of equations 
Re A = Re B — and thus takes the integral form, 

Im (f> w(z) dz = 0. (4.6) 

Ja,b 
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Equations (4.6), known as the Boutroux equations, are the substantial ingredient 
of our proof below of the solvability of the RH problem. 

Denote D = z\z\ the module corresponding to our elliptic curve. We ob- 
serve that ( |4.6D determines D as the diffcrentiable function D(ip) of the angle 
parameter ip £ (^,7r). Indeed, the l.h.s. of ( |4.6D yields the complete set of the 
independent first integrals of the first order nonlinear ODE 



dD _ 1 

dip 41m (LO a u)b) 



(Cu a (e^n a - e'^n a ) - u b {e^n b - e-^0 6 )), (4.7) 



where 



dz 

a.b W{z) ' 



z 2 dz 
b w ( z ) ' 



Whenever the elliptic curve w(z) does not degenerate, the r. h.s. of (4.7) is a 
smooth function of ip, and, given the initial point, equation (4.7) is uniquely 
solvable. To find the initial point, observe that equations (4J3) admit degener- 
ation Z\ = Z3 fo r ip = it and z\ = for ip = At these boundary points, 
the r.h.s. of (4.7) remains continuous which ensures its solvability. The classical 
uniqueness theorem for ODEs breaks down at these boundary points, however 
the only one of the solutions of the differential equation ( p~7|) passing through 
the initial boundary point satisfies the integral equations ( |4.6| ). 

Define the Stokes line for the linear equation ( |2.6| ) emanating from the point 
Zk, k — 1,2,3,4, as the set of the points z satisfying the equation 



Im 



w{z) dz = 0. 



(4.8) 



Define the Stokes graph as the union of all the Stokes lines. The conditions (4.6) 
mean that the Stokes graph is connected. On the other hand, conditions (4.6) 
allow us to choose the arcs of the Stokes lines as the segments [23, 21], [zi, — zi] 
and [-21,-2:3]. 

For the following, we need to know the jump conditions of the function g(z). 
These are described by the equations 

1 



z £ {z a ,zi): g+ + g- = 



z £ {zi,-z x ) 



9+ 



-A: 



z £ (-zi,-z 3 ): g + +g_ = --B. 



(4.9) 



5 The model Baker- Akhiezer RH problem. 

We remind that the reason why we replaced the original RH problem by the 
one presented in Figure 10 is that 

e -taW<ra Se tg{*)*» j ; (5J) 
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as \t\ — > oo for all z belonging to the infinite branches of the jump contour 7 
and z ^ ±Zfc. Indeed, equation (5.1) means that when we normalize the RH 
problem, i.e. when we pass from the function '5 to the function 

Y(z) = ^{z)e t9 ^ a \ 

so that Y(z) — ► / as z — > 00, the Y"-jump matrices will exponentially fast 
approach the identity matrix as |i| — > 00 everywhere on the infinite branches of 
7. Therefore, one can expect that the main contribution into the asymptotics of 
the solution of the RH problem will be made by the truncated RH problem posed 
on the sum of the segments (see Figure 6) (2:3, z\) U {z\, — zi) U (—Zx, —23). In 
other words, we arrive to the model problem consisting of finding the piecewise 
holomorphic function ^ BA ^{z) with the following properties: 

1) near infinity 

¥ BA \z) = (I + 0{z- 1 ))e' tea3 , z^oc, (5.2) 
where 6 = i(fz 3 + e tip z); 

2) on the sum of the oriented contours (23, Zi)\J (z%, — Zi)U (— z\, — Z3), the jump 
conditions hold 



z G (-z u -z 3 ): * ( + A \z) = y ( _ BA) (z)S A , S A = ( ^ •» ) . (.-,.:■{) 



zefazx): *V ; (*) = *- (*)Sfl, Sb=I k _ ± "J, (5.4) 

*e(*i,-*i): * { + BA) (z) = ^ BA) (z)Sc, S C = (j=*Z !_ SlS3 ); (5-5) 

3) the function ^ BA \z) has no more than L2-integrable singularities at the 
points ±zi, ±23. 

We use the superscript (BA) since we are going to construct ^ BA >[z) in 
terms of the Baker- Akhiezer functions of the Riemann surface T. On this Ric- 
mann surface, there exists the only one, up to an arbitrary constant multiplier, 
holomorphic differential with the periods uj a ,b- Introduce the canonical 
Abelian (elliptic) integral of the first kind U(z), 

U(z) = - f -^y (5.6) 



and the Abelian integral h(z) which differs from g(z) (4.4) in some integral of 
the first kind, 

h(z)=g(z) + cU(z), (5.7) 

where the constant c will be determined later. Since U(z) = 0(z~ 1 ) as z — > 00, 
the integral h(z) has the same canonical asymptotics as the integral g(z), 

h(z) =i(^-z 3 + e lv z) +0(z- 1 ). (5.8) 



20 



Denote the corresponding complete integrals by the letters A' and B'\ 



A' = U (f w(z) dz + — I = A + c, 

J a U a J a w(z) 

B' = Ai J> w(z) dz+ — <L -^rx = B + ct, t=— . (5.9) 

Jb U a Jb W(Z) UJ a 

o 

Let us introduce the "intermediate" i]/-function, 

^{z) = ¥ BA) {z)e toh{z),J \ (5.10) 

where to is the second reserved parameter to be defined later. The asymptotics 
of this function at infinity, 

^{z) = {l + O{z- 1 ))e- {t - t0 ^' y \ z^oo, (5.11) 

differs from the canonical in the difference t — to which replaces the original t. 
Its jump properties are similar to those of the function vI/( Bj4 ) except for they 
are described using the new parameters to, A' and B': 

z e {- Zl ,-z 3 ): 9 + (z)=&..(z)Ga(z), (5.12) 

ze(z 3 ,zi): &+(z) =&-(z)Gb{z), (5.13) 

ze(zi,-zi): & + (z) = $-(z)G c (z), (5.14) 

where 

A ( -s?e*° B '/*\ , _( S3 e^B'/2 
- \s 3 e^ B '/ 2 ) ' &B ~ {-s-^ B '/ 2 

(1 - SiSz)-^-* ^ 



Gn = 



(l-«i* s )e*° 



t A' 



Taking into account the conditions (3.11), we may choose the reserved pa- 
rameters c and to in such a way that 



*3 



= e tnS '/ 2 , l-sisa = e~ toA '. (5.15) 



The system (5.15) yields immediately the parameter values 

t = -2— In s 3 - —ln(l - sxs 3 ) (mod - ) (mod - ), 

c n A, B, , . . A-rtiA^ , , 2iriB . , „, 

— to = 2-lns 3 + -ln (1 - «i« 3 ) mod—— mod——, 5.16 
u a A A A A 



where 



A = Auj h - Bu a = 4i(j a uJb - JbU a ) = ^~( z i 



2 J, £) = ^ e i<p 
3) 3 
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due to the Riemann bilinear identity. 

, . o 

Under the conditions (5.15), the RH problem for the function simplifies 
dramatically: 

1) near infinity, it behaves like the exponent 

$(js) - (7 + 0(z- 1 ))e~(*- t °) eCT3 , z^oo; (5.17) 

2) it jumps on the intervals (—21, —23) and (23, 21): 

ze{-z 1 ,-z 3 ): y + (z)=$-(z)( 1 (5.18) 
2 G (23,21): | + ( 2 )=i(z)|_ j (5.19) 

3) it has no more than L2-integrable singularities at ±21, ±23. 

Apart from the countable set of values of t described below, this problem 
can be solved explicitly, cf . Jt| . Let us introduce the function 

a- l -( ^ + -^~^)\ (520) 

where 

/ r 2 - 7 2 \ 1/4 

is defined on the upper sheet of the Riemann surface T (see Figure 10) in such 
a way that 

(3 — ► 1 as 2 — > 00, 

so that 

/i — ► /. (5.22) 

z — >oo 

The discontinuity of the scalar function (3 is described by the relations 

2 G (-21,-23): P+ = -i0-, 
z G (23, 21): /3+ = 

and therefore, for the function /x, the jump conditions hold: 
2 G (-21,-23): /i + =/i_^ 1 

2 G (23,21): /i + = /i_^_ 1 (5.23) 

o 

The function \f» is constructed in terms of the Baker- Akhiezer functions as fol- 
lows: 

k*)=( fXll l BA M12 |!" A | ^ 
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where 



and 



@ _ - (t -to)(n W -^) e(u(z) + v(t-t Q ) + ^)e^) 

^ BA[) e(u(z) + ^)e(v(t - 1 ) + 1) ' 

*^- e e(t/(z) + i)e(^-i ) + i)' (5 5) 



Q(z) = <?(*)- AE/(z), U{z) = —( Z ^- r=^, (5.26) 

V — ^ fx 2 + z 2 ) - e iv 

27riw a 3w a 3 3uj a ' 



and 6(z) means the Riemann theta-function of T, Q(z) = e 



Tzin T+2-Kinz 



(k) (fc) 

It is worth to note that all vf ba(z), ^*ba( z ) are single valued on C * cut 
along [z 3 ,zi] U [—21,-23], see Figure 10, because § a dVL — 0. Also, considered 

(fc) (fe) (fc) (k) 

on the Riemann surface, $si(z) = vf ba{P) and ^* BA (z) = vfi Bj 4(P*). 

Now, let us verify that the function (5.24) yield the solution of the RH 

problem ( 5.17 )— ( 5 ■ 18| ) . Note first of all that ^(2) has no pole on the upper 
sheet of T. Indeed, Q(U(z) + 5) ^ Vz £ C*. Taking into account the equation 
6(i±i) =0, the only zero of 6(C/(z) + i±^) lies at 00+. Thus 6 (U{z) + ^f) = 
0(z~ 1 ) while both fii 2 and [121 vanish as 0{z~ 2 ). 

o 

Next, to check the asymptotics of ^f(z) near infinity on the upper sheet of 
r, we find 

U(z) = — + 0{z- 2 ), z^oo+, 



LO n Z 



and therefore 



' 22 

Similarly, using the asymptotics 



(*) 11 = e - i ( t - t 0)(^ 3 + e< ^)(l + O(2- 1 )), 

(&) = ^-t o) {^+^ z ) ^ + o{ z - 1 )) . 



Vi2 = -i Z \Jt + {z-% ^ 1=l i-^l+0(z- 
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we obtain 

(*) ai = e -(*-*o)(!* 3 +^*) . {z- x ), (#) 12 = e ^(*-*o)(^ 3 +e-,) . 



and the asymptotics (5.17) holds true. To check the jump conditions (5.18), 
( |5.1E ), it is enough, due to (5.23), to show that 

(fc) (k) 

on the segments (z 3 , Z\) and (— z%, —z 3 ). Let us consider the first of them. One 
can see that on the segment (z 3 , Z\), 

U + + U- = T -, Q + + Q_ = ^(B-A T ) = -^- = -mV. 

Therefore 

w+ c - f t- fa ^ e(M*) + y(t- to) + l)e(|) _ 

= (t -to)(n-(z) +ff «v) + I + - to) + |)Q(|) _ 

e(-[/_( 2 ) + 5 + |)e(F(t-t ) + |) 

_ e ( t - t0 )(n_ (z )-,^) e(^(z) + l±i-y(t-t ))6(l) _ m 

9([/l(*) + ±^)G(F(i - to) + \) * ba[z) ~> 

where we have used the properties 0(— z) = O(z), 9(z + 1) = O(z), 0(z ± t) = 
e _7rirT27rM 9(z). The other conditions can be checked similarly. 

Thus we obtain the following representation for the function ^^ BA ^: 

( (1) 

= | e -*0^3 = W*BA ] e -toha, = (5 _ 28) 



e (* J 4-t A')(7cr3 e -i g cr3 



where 

q)_ e(C/(z) + T/(t-t ) + i)9(l) m_e(U(z)-V(t-t ) + ^)e^) 



e(c/(z) + i)e(F(t - to) + 1) ' e(f/(z) + ^)e{v(t - 1 ) + ±) ' 
m_ Q{u( z ) + v{t-t ) + i¥)Q{±) ™_ Q{u( z )-y(t-t Q ) + \)e{\) 

G(U(z) + ±¥)Q(V(t-t ) + ±)' B(U(z) + %)9(V(t-to) + %)' 

Note finally, that the determinant of ^ BA \z) might have singularities at 
the branch points ±zi, ±23 of order 0((z — z/j) -1 / 2 ). However, because of 
the asymptotics det^ ( - BA ^ (z) = 1 +0( z X J as z — ► 00 and absence of jump, 
the determinant is entire function of z which, due to the Liouville theorem, is 
constant, det^ (BA) (z) = 1. 
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6 Local RH problems near the branch points 



6.1 RH problem solvable by the Airy functions 

The function &( BA '(z) has branch points at z — Zk, k — 1, 2, 3, 4, and does not 
approach the solution ^(z) of the problem ( 3.28| )-(3.35). In this section, we 
construct such an approximation using the model problem related to the Airy 
functions. Consider the Wronski matrices of the independent solutions yi(Q of 
the Airy equation y" = (,y: 



*(0 



2/i (0 2/ 2 (C) 
2/i(C) 2/ 2 (C) 



Using asymptotic formulae presented in ref. 0, one can introduce the matrix 
functions $fc(C) with the standard asymptotic behavior as £ — > oo, 

**(C)=r*(} ^(i + Or^^ef^, C 3/2 -oo, (6.1) 
Ce^ = {(eC: argCG (^(fc-|) + e , ^(fc + Ij-c)}, e >0. 
These functions are related to each other by the Stokes multipliers 

$ fc+1 (C) = <i> fc (C)G fc , Ga*-i = (j ~*V G 2fc =^. JV (6.2) 

This collection of the matrix functions allows us to construct an exact solution 
of the following RH problem (see Figure 12): 

1) near infinity 

* = r*(j \)(l + 0{C m ))e^ l2 °\ C 3/2 -oo, (6.3) 

where £ 3 / 2 is defined on the plane cut along the ray arg£ = — 7r/3; 

2) on the anti-Stokes rays arg£ = 4^fc, k — 0, 1,2, oriented to infinity, the 
jump conditions hold, 

$+ = $_G fc , argC=yfc, k = 0,1, 2; (6.4) 

3) on the ray arg£ = — 7r/3 oriented to infinity, the jump condition holds, 

$+ = $_Goo, G oc =ia 1 =(. *Y (6.5) 
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Figure 12: The RH problem graph for the collection of Wronski matrices of the 
Airy functions. 

6.2 Parametrix in a neighborhood of the branch point z% 

Let us change the variable in accord with 

C(z) = (-^9s(z)) V \ 93 (z) = 4iJ w(z)dz = g(z)--B. (6.6) 

This change is defined on the complex z-plane cut along the sum of intervals 
[zi, —zi]\J[—Zi, — z 3 ]U [— z 3 , oo) and mapping the domain T> 3 of z-plane bounded 
by the Stokes line chain shown in Figure 13 onto the complex £-plane cut along 
the ray arg£ = — tt/3. Now, the function 

$(z) = $( C (z)) (6.7) 

which behaves at infinity as follows, 

*(*)=(-§<73(*)r ¥ (} -l)( / + °(t^))) e " tflsWOS ' ^sW-oo, 

(6.8) 

jumps on the graph shown in Figure 13. 

The rays arg£ = 2-7rfc/3, k = 0,1,2, are transformed into the anti-Stokes 

A (3) „ (3) ^ (3) 

lines 75 , 7g and j\ , respectively, see Figure 4, moreover the jump matrices 
are preserved there. The cut arg£ = — n/3 turns out the segment (z 3 , z\) which 
splits at the point z\ into the boundary of the ^-domain. 
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Therefore, the function 



$( 3 )( z ) = $(z)e^ CT3 s~ 



(6.9) 



has the j umps o n th e graph shown on Figure 12 which is described by the 
matrices ( ^30| ), 



(3). c (3) 



1 

i 1 

S3 

1 

i 1 

S3 



z£lf ] : S 6 



z € (z 3 , zi): Si 



1 -s 3 
1 



_i_ 

S3 



and has the asymptotics near infinity 



(6.10) 



*53(^) 



)) 



(6.11) 
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It is important, that the jumps of the matrix functions \t (z) and <I>( 3 ) (z) coincide 
with each other on the contours 7ej 3 \ and (23, Z\), so that they differ 
from each other in a left matrix multiplier (z) holomorphic in some bounded 
neighborhood of the node point z = z%: 

*(z) = M {3 \z)^(z). 

Besides the Airy fundamental matrix set let us introduce its reduced 

asymptotic form 

*«(c) = r*(; _ 1 1 )e*^, 

which solves the RH problem on the ray arg£ = — 7r/3 oriented from zero to 
infinity, 

$a S+ =$a S _Goo, G ao =ia 1 = (, *V (6.12) 



Using the changes ( |6.q ), ( J6.9D , the matrix 



jumps on the Stokes segment (23, zi), and the jump matrix is (3.34) 



s B = 



S3 



z e (23, zi). 



Since the jumps of the matrix functions ^^ BA \z) and $o» (2) coincide with 
each other on the segment (23, z\), the functions differ from each other in a left 
matrix multiplier N^(z) holomorphic in some bounded neighborhood of the 
point z = 23, 

¥ BA \z) = N (3 \z)^Xz)- (6.13) 



Let us define 



(6.14) 



6.3 Parametrix in a neighborhood of the branch point Z\ 

To construct the parametrix near z = z%, let us introduce the function ^(^(O, 



$(D( C) 



*o(CK 

*i(0 



. 1-S1S3 



r for argC6(-f,0), 
"(i^)"* fOT ar gC e(0,f), 
'(i^Sr)"* for argCe(f,7r), 



(6.15) 



$2(C)e^ fCT3 (3 3 (l-3i3 3 ))' 
[$3(C)e-^ CT3 (ss(l-*i«3))' 



for argC € (tt, 4p) 



3 

for ar gCe(f,f) 
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As it is easy to see, the introduced matrix function satisfies the following jump 
conditions on the rays oriented to infinity, 



4£ 


(0 




c(l) 

°1 


*£ 


(0 




c(l) 




(0 




S c = 




(0 


= $ w (C)s«, 


c(l) 
°5 




(0 


= <*> w (C)s B \ 





S3 

1 



I-S1S3 

f 



argC = 0, (6.16) 



argC = y, 



l— si«3 



1 - S1S3 

1 



S3(l-SlS 3 ) 

-S3 



argC = 7T 
argC = 



47T 



argC 



7T 

3' 



Introduce the change of variables 

C(z) = ( 

where 



3i 



,2/3 



2 Pi (2))' . 
51 (z) =4i w(z)dz 



(6.17) 



is defined on z-plane cut along (00, 2:3) U (23, Z\) U (00, —z\) U (—21, —23). This 
change is holomorphic near the point z\, and hence it is defined on the 2-plane 
cut along the sum of intervals (23,00) U (00, —21) U (—21,-23) and maps the 
domain T>\ of 2-plane with the boundary on the Stokes lines shown on Figure 14, 
onto the complex (- plane cut along the rays arg £ = — w/3, arg(£ — ( 2 ) — tt, with 
some (2, arg ^2 = t. Thus, in the area to the right of the cut (23, 21) U (21, —21), 
see Figure 13, or, cquivalcntly, for arg(— ^gi(z)) € (— § , ^ L ), 



ffi(s)=fl(s)-^-^4 



(6.18) 



and in the area to the left of the cut (23, 21) U (zi, —21) i.e. for 2 satisfying the 
condition arg(-f 9l (z)) e (4f , 4f), 



5l (2)= 5 (2)-^+-A 



The function 



$(!)( z ) = $«(C(2)) 

with the asymptotics at infinity 



(6.19) 
(6.20) 



3f 



[t 3 e i(tB-t B')a 3 e i(tA-t A')a 3 
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tgi(z) -> oo, arg(-%(z)) e (-|, y), 

xe -tffo-3g-if o-3 e 3(tS-t S')°'3 e -5(tA-toA')cr3 (6.21) 

igi(z) -> oo, axg(-tgi(z)) e (y , -y)> 

jumps on the graph shown on Figure 14. Under the change ( |6.17 ), the interval 
(z 3 , zi) corresponds to a segment of the ray arg£ = — ^ (if the orientation of the 

ray is from infinity to the origin), the curve j[ turns into the ray arg£ = 0, 
the curve 72 transforms into the ray arg£ = the Stokes line {z\,— z\) 
corresponds to the ray arg£ = 7r and the curve corresponds to the ray 
argC = 
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Since the discontinuities of the functions ^(z) and ^^\z) are the same on 
the contours 7^ and the Stokes segments (z 3 , z\) and (z\, — z\), two these func- 
tions coincide with each other up to a left matrix multiplier M W (z) holomorphic 
in some bounded neighborhood of the point z\\ 

Besides the model function $W(z), we introduce its reduced asymptotic 
form $1V (z) as follows: 

*i^)=(-f^))~"(j \ 

xe -tg^3 e -i^cr 3 e ^(tB-t B' )a 3 e ^(tA-t A')a 3 ^ 

tgi(z) -> 00, arg(-t 9 i(z)) e (~, y), 

^w=(-ffliw)-*(; \ 

xe -tgtT3 e -i^C3 e i(tB-t B')(T3 e -^(tA-toA')cr3 ^ (6.22) 

tgi(z) -» oo, arg(- 9 i(z)) £ (y , y), 

This reduced function is continuous across the anti-Stokes curves 7^, and is 
characterized by the jump conditions 

*£ + (z) = ^_(z)Sc, S c =(^ l- ai 8 3 )> * e (^-*b>> 

Note again, that the jumps of the functions ^^ BA \z) and <&£){z) coincide on 
the segments (23, z\) and (21,-21) and hence the functions coincide with each 
other up to a left multiplier N^(z) holomorphic in some bounded neighborhood 
of the point z\ : 

^ BA \z) = N^\z)^}{z). (6.23) 

Define 

& 1 \z) = N( 1 \z)& 1 \z), NW(z) = & BA \z)($W{z)y 1 . (6.24) 

The functions solving the model RH problem associated with the branch 
points —z\ and —z 3 can be constructed in the very same way, however they can 
be easily defined by the use of the symmetry relations: 

¥ 2 \z) = a 2 ¥ 1 \-z)a 2 , tfW(*)= o- 2 tt (3) (-*)<r 2 . (6-25) 
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7 Asymptotic solution of the main RH problem 



In this section, we will construct the approximate solution ^(z) of the RH 
problem fl3.28|)-([n4|) and prove the solvability of the latter. To this end, let us 
introduce four not overlapping disks Bk, k — 1,2,3,4, centered at the branch 
points Zfc, circles Ck = dB^, and define the piecewise holomorphic function 

, ( ) _ r*W(z), zGBk , k= 1,2,3,4, , v 

The function "J solves the RH problem on the contour depicted on Figure 15: 



C 



3_ 



C, 



-z 



Figure 15: The model Riemann-Hilbert graph. 



1) near infinity, due to (5.2), 

= ¥ BA Hz) = (I + 0(z- 1 ))e- te ' 73 1 z 



oo, 



(7.2) 



where = i(f z 3 + e^z); 

ffel (k) 

2) on the parts 7, H Bk of anti-Stokes lines 7,5 enclosed in the interior Bk 
of the circle Ck, and on the segments (23, z\), (zi, —z\), (—21, —23), the jump 
conditions hold 

= *_(z)S (7.3) 



with the jump matrices S defined in ( 3.30 )-( 3.34 ); 

3) on the circles Cfc, fc = 1,2,3,4, oriented counter-clockwise, the jump 
conditions hold 

*+(*) = (7-4) 
where, in accord with (|6T3|) , (fTlJ), |[23|), (^2^), and d^H), 

iJ(z) = = (^]{z))- 1 ^ k \z), z e C k . 

To describe the difference between the model problem solution \P(z) and the 
solution ^f(z) of the original problem, let us construct the "ratio" 



Z(z) = y{z)y-\z). 



(7.5) 
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The correction function Z(z) satisfies the RH problem on the jump contour 
shown on Figure 16: 




Figure 16: The RH graph for the correction function Z(z). 



1) Z(z) is normalized at infinity 

Z — ► I (7.6) 

2) on the parts of the anti-Stokes lines outside the circles Ck and on the circles 
themselves, the following jump conditions hold true: 

Z + {z) = Z.{z)Gf\z), ze^n(C\B k ), (7.7) 

where, on the anti-Stokes tails, 

G « (z) = $(^)( z )s(»)(f(^)( z) )- 1 i 

and, on the circles Ck, k = 1, 2, 3, 4, 

Z+(z) = Z_(z)G k (z), zeC k , (7.8) 

where 

G k {z) = ¥ BA \z)(¥ k \z)y l = 
33 



= 7VW(z)$( fe s )(z)($W(z))- 1 (^ fe )(z))- 1 . 



As it is easy to see, for t lying apart from the points (5.27), the jump matrices 
on the anti-Stokes lines have the asymptotics 

G ( s k) (z) = I + ( e - 2 l Rc (*9( z ))l) =/ + ( e - 2 l Ro (*9f=(z))l) ; 



where we have used the Boutroux equations (4.6). For t lying apart from the 
points (5.27), using the definition of (z) and the boundedness of the matrices 



N(k)(z) and their inverse, we find the estimate 

G k (z)=I+(t- 1 g^(z)). 



Now, we are going to prove the solvability of the RH problem (7.6)-(|7.8|). 

First note that the jump matrices are continuous on the RH problem graph. 
Consider the connected component of the graph associated with the point Z3 
(see Figure 17). 




Figure 17: The connected component of the RH problem graph associated to 
the point Z3. 



Seeking for convenience, we orient the anti-Stokes lines 7^ in the direction 
from infinity to Z3. Denote provisionally the jump matrices for the function 
Z(z) as R(z), i.e. Z + (z) = Z-(z)R(z), moreover: 

(3) 

on the part C 25 of the circle C3 between the points of its intersection with the 
Stokes line (z 3 , z{) and anti-Stokes line the jump matrix R(z) — Rq(z) is 
as follows: 

i? (z) = ^ (3) (z)^(^^ (7.9) 

~ (3") 

on the anti-Stokes line 75 , the jump matrix R(z) ~ Ri(z) is given by 

R x {z) = ^ {BA) {z){Sf ) )- 1 {^ BA \z))- 1 = I + 0{e 2t9(z) )- (7.10) 
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(3) 

on the part Cg 6 ; of the circle C3 between the points of its intersection with the 
anti-Stokes lines 7^ and 7g , the jump matrix R(z) — Riiz) is as follows: 

R 2 (z) = N^{z)^}{z){$f \z))-\N^(z))- 1 = 

= M 3 \z)^(z)(s^)-\^\z)yHN^(z)r 1 = 

= Ri(z)Rq(z) = I + 0(t^ 1 g^ 1 (z)); (7.11) 
on the anti-Stokes line 7g , the jump matrix R(z) = R^{z) is given by 



(3) 

on the part Cg-/ of the circle C3 between the points of its intersection with the 



R 3 (z) = ^ BA ^(z)(S 6 y 1 (^ BA) (z))- 1 = 1 + 0( e - 2t9(2) ); (7.12) 
)f the circle ( 

anti-Stokes lines 7g and 7J , the jump matrix R(z) = R^{z) is as follows: 

R,(z) = N^iz^Jiz^iz^iN^iz))- 1 = 

= M 3 )(z)$( 3 s )(z)(S 6 )- 1 (si 3 V 1 K 3) (z))- 1 (iV( 3 )(z))- 1 = 

= R^R^Roiz) = I + Oit-ig-^z)); (7.13) 

on the anti-Stokes line 7 J , the jump matrix R(z) — R§(z) is given by 

R 5 (z) = ^^ BA \z)(S < f ) y 1 {^^ BA) {z))- 1 = 1 + 0(e- 2tg(z) ); (7.14) 

on the part C{ 2 of the circle C3 between the points of its intersection with the 
anti-Stokes line 7 J and the Stokes line (23, zi), the jump matrix R(z) = Rq{z) 
is as follows: 

R,{z)=N^{z)^{z)^f\z))-\N^{z))- 1 = 

= N^\z)^(z)(sfYHS e )-\strH^\ Z )rHN^Hz))- 1 = 

= R 5 (z)R 3 (z)R 1 {z)R (z) = 1 + Oit^g^izj). (7.15) 

The relations 

i?2 = RiRq: Ra = R3R2, Re — R5R4 (7-16) 
above provide us with the continuity of the RH problem at the node points 

(3) A — I 

C3 n % (see |27| ) . The fact is obvious because of the possibility of analytical 

( 31 

continuation of the jump matrices R(z) from the contours C3, 7£__and the 
possibility to split the graph at the node points in accord with ( 7.16 ) and to 
transform it into the union of the circle C3 and of three smooth curves with the 
common endpoint P = C3 n (23, z\). As to the point P itself, the continuity 
of the RH problem follows from the equation Ro(z) = Re(z). Indeed, at this 
point, 

Ro(z) = ^(z^iz^iN^iz))- 1 = 

= ^ BA \z)S B (^\z)S B )- + \N^(z))^ = 

= ¥ BA \z)(^\z))-\N^(z))-i = R 6 (z), (7.17) 
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where plus and minus mean the left and right limits on the segment (2:3, z\). In 
the very similar manner, one can check the continuity of the RH problem on 
the other connected components of the graph. 

All the jump matrices are uniformly close to identity. More exactly, for t 



lying apart from the points (5.27), the estimates hold true 



( „ „-2|Re(tg(2))| _ c ~(*0 

w*)-n<\ Cl ?Ju \-i z % ls r > (7.18) 

l ci\t\ L \gk{z)\ L , ze C k , 

where the concrete value of the constant c\ is not important for us. The esti- 



mates (7.18) imply that the differences R(z) — I are square integrable on the 



RH problem graph E = U k C k U k , s 7? } , i.e. R{z) - I € L 2 (£) 



The RH problem (7.6)-<XLH), i.e 



Z — > /, 

z—>oc 

Z + (z) = Z-(z)R(z), zeS, (7.19) 
is equivalent to the system of singular integral equations for the function = 

P(z) = I + I p(0(R(0 - I)t^~, (7-20) 

where z_ denotes the right limit of z on the contour of integration. In the 
symbolic operator form, the system reads 

p = I + C-[Mp], (7.21) 

where C_ is the Cauchy operator, and M is the operator of the right multipli- 
cation in the matrix i?(£) — I £ L 2 (Yj). Thus M acts in £ 2 (£) and satisfies the 
estimate 

||M|U 2(S) < c 2 \t-\ (7.22) 

where the precise value of the constant C2 is not important for us. Because the 
Cauchy operator is bounded in £ 2 (S), see |24|, 1 1|, we obtain the estimate 

||C_M|| L 2 (E) <c 3 |t- 1 |. (7.23) 



Therefore, for any large enough t lying apart from the points (5.27), the operator 
C-M is contracting, equation for the function \ = p — /, 

X = C_[M]+C_[M X ], (7.24) 

is solvable in the space L 2 (T,) by iterations, and the estimate holds true, 

\\p-I\\ L 2 (j:) <c\t\-\ (7.25) 

The solution of the RH problem ( 7.19| ) for the function Z(z) is given by the 
Cauchy integral 



Z(z)=I+^ I p(C)(R(C) - I)- 
2m 7 S ( 



3G 



which implies the asymptotic relation 

*(z) = (/ + 0((l + |z|)-4- 1 ))*(z) 
in any closed subdomain of C \E. 



(7.26) 



8 Asymptotics of the Painleve function 

Now we are ready to compute the asymptotics of the Painleve function. Using 
the expansion of ^(A) near infinity, 

*(A) = (J + \(~Da 3 + \ua x ) + 0(A)) e -(t* 3 +* = 



where u is the Painleve function (cf. (2.7)), i.e. the solution of (1.1), D = 



= (7+ 0(— ))*(*) = (/ + Oit-'z- 1 )) y(z)e- tah{z)a \ (8.2) 

to find the asymptotics of u up to terms of order t , it is enough to expand 
the Baker- Akhiczcr function <f Bj4 or Sj 4 near infinity: 

2tV3 - A Ze 9 (t/(z) + ^)e(V(t — to) + \) ] 

= iOztf _ 2)e ^- t0) y e(y(t - to) + ^)8(|) 

4 V1 3; e(y(t_t ) + i)e/(i±i) v ^ v ^ 



E quat ions (4.6), ( 5.16 ), and ( |8.3| ) are transformed to the statement of The- 
orem [2.l| via routine manipulations with the elliptic functions and integrals. 
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